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ABSTRACT 


A digital computer simulation program using Digital Simulation 
Language is produced to study the performance of a hydrofoil in calm 
water. Various automatic control systems are studied with the model 


constrained to the pitch-heave-surge mode of operation. 
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I. INTRODUCTION 


A, COORDINATE SYSTEMS 

Three coordinate aa will be used in this study of hydrofoil 
craft dynamics. They are: the body axes, the water axes and the earth 
axes. Each is a right hand orthogonal system. 

The body axes coordinate system (Xo Yo a Powe in figure 1-1 
has its Origin at the center of gravity of the craft and is fixed relative to 
the craft. This system is initially coincident with the origin of the earth 


§. : 
axes system, however, it moves with the craft as time progresses. X is 








Poeitive P 
(Rolling t> St) 


2 cy Ponitivo Q 
<=. 
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(Yewing to Stbd.) 


‘ Positive W 


Orientation of the Body Axes with Respect to the 
Craft and Directions of Positive Velocities 
Figure l]-] 


positive forward, Y is positive to starboard and Z is positive 
downward. 

The earth axes coordinate system XK, Nae, 2) is fixed relative to 
the earths surface. The origin may be chosen at any convenient location, 


provided the body system origin is initially at the same point. The x, 


li 





axis lies in the horizontal plane; it is initially coincident with, and 
positive in the same direction as, the body X axis. The <— axis lies in 
the horizontal plane and is positive to starboard when the Observer is 
facing the positive Xp direction. The Za, axis is normal to the horizontal 
and is positive downward. Water surface motion and crafts' motion 
relative to calm water are described in this coordinate system. 

The water axes coordinate system is aligned with the relative 
velocity vector and resolves into ‘lift’, 'drag' , and ‘side force' directions. 


All hydrodynamic model test data is taken and presented as plots in terms 


of this coordinate system. 
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Definition of Positive Directions for Euler Angles 


Figure ]-2 
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Detailed derivations of the transformation matracies between coordi- 
nate systems are Outlined in Ref. 1 , however, only the results will be 
considered here. In general, the body axes are displaced from the earth 
axes by the Euler angles, OD ,©C , wamid ve . Figure 1-2 defines the posi- 
tive directions for the Euler angles and established the means for trans- 
formations between earth and body axes. 

The hydrodynamic force and moment data are Obtained in a water 
axis system which is always oriented with respect to the relative water 
velocity. These data must be transformed into body axes so that their 
effects can be included in the equations of craft motions. In general, the 
water axis Orientation is always changing with respect to the body axes. 
This necessitates a transformation to resolve lift, drag, and side force 
quantities into body axes. The orthogonal water axis coordinate system is 
defined in figure 1-3. The angles of rotation between water and body axes 


are & and 8 and will be defined later. 


LIFT 


SIDE SLIP 


DRAG 


Figure 1-3 Definition of Water Axes 
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B. EQUATIONS OF MOTION 

The development of the six equations of motion for a hydrofoil 
follow the same classical lines as those for a displacement type hull. 
This development is clearly outlined in Ref land4_ and therefore 


will not be done here. The six equations for the hydrofoil are: 


Translation 


IL= mk -QW + RV (1-1) 
| t VemR+PW-RU | (1-2) 


: | 
\W = 2tQU- PV (EC) 


Rotation 


P=,,[L-QR(Ing-Ty)+(R+0P) Zea) WS 
Q= Z, [= RPI = ~(p-R) I, | Gee) 


R=TAIN-PQ(IyTy)-(or-P)Ze] ave 


Comparison of actual and simulated data show that the translation 
terms RV, PV, and PW may be neglected. Similarly, the product of 
inertia quantity, Leg! is less than 10% of ix and 4% less than Loe" 


Consequently, when considered in conjunction with the magnitudes of 
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(R+PQ), Woes and (QR-P), the terms containing Lg may be neglected 
fOnsthisucnati- 

Given the equations 1-1] through 1-6 for linear and rotational 
accelerations in body axes, the velocity terms in body axes are obtained 
by merely integrating each equation one time. These outputs must then 
be transformed to earth axes to be of useful form. The transformation of 


the linear equations become: 


Ue = UCOS 6 Cos ¥ +V(Cos SING SIND - SIN Y SIN @) 
+Wcos ¥ Cos 6 Cos @+SINY SiN g) (Ie7) 


Ve = USin ¥ CoS O+V(COSY¥ Cos P+ SIN SIN © SIN O) 


(=3) 
+W(sin YSIN@ GOs0 = Gael Sik ¢) 


We = ~USIN@ +VCOSOSING + WLoS 6 Cos ¢ (1-9) 


Only the crafts velocity in the downward direction is required to 
compute the crafts height above the water surface. Therefore, the verti- 


cal position in earth axes is given by: 


Ze> |W, dt (1-10) 


The crafts position in earth axes coordinates is found from: 


Mee [Up dt Yen | Vedt (1-11) 


The Euler angles, a) ,8, and¥Y ,. are also of interest and can be 


found by integration of the Euler rates cd ee anc ye : The rates in 


turn must be derived from the body angular rates. The Euler rates are not 
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easily obtained because they occur about axes which are not orthogonal. 
This fact can be appreciated by referring to figure 1-2 and recalling how 
the angles were defined. ai is a rotation about the Zo axis, QO isa 
rotation about the XY) axis, and Dp is a rotation about the X axis. The 


results in the following set of equations: 


g = P+ YSIN 6 (3-12) 
6 =Qlosgd-RSiNG (1-13) 
Y = (Qin d+ RCosd)CosO +($-P)Sin © (1-14) 


Studies of this craft show that even under the most drastic abnormal 
situations, pitch angles are not expected to reach 10° and under normal 
operating conditions, roll angles will not exceed 3°. Therefore, the 
following assumptions are justified: 


Sin@ = © (1-15) 
Cos@= 1 


The Euler angles are obtained by integration of equations 1-12, 


13, 14. 


©. CALCULATION OF VELOCITY COMPONENTS 

The variOus foil velocity components must be determined at this . 
point in order to calculate the angles of attack and side slip. The total 
velocity components of each foil in body axes can be expressed in terms 
of the craft linear and angular velocities by considering the craft geometry 


as represented in figure 1-4 and Table l. 
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The resultant foil relative velocity equations become: 


Center foil 


U= Ut _ Q (1-16) 

Ve=V— Lovgeneer L een (1-17) 

Wet Wo Lycee Q (suey 
Port foil | 

Up= Ut Lao Q-Lype R (1-19) 

Vp=V- Lape Pth ype R (1-20) 

WpetW-LyprpP—LyppQ ae 
Start Oard foil 

Ug Ut LageQ-Lyoe R (1-22) 

Vee Vo Leet t Lye R (1-23) 

Ws =W+ Lyse P- Lge @ (1-24) 
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Pigure 1-4 Definition of Symbology 


© = point of application of foil forces 

© = p0dint of application of strut forces 
Note: 

xcr = ‘xcs 


Lypp = Lyps = Lysp = *xss = ‘koe 
“YSF 7 Lyss ~ ~lypr = “typ 


Lopr = “sp > ‘ome 


Lops 2 Loss are variables and in general will not be equal. 
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Mid foil 


Ue iO (1-25) 
Vie = oa Lewes f? tloyme R (1-26) 
Wy = Wo Lane @ (1-27) 


These foil velocities are in body axis terms and can be used directly 
to calculate the angles of attack and angles of side slip. 

The angle of attack of the foil is shown in figure 1-5a and defined 
as the angle whose tangent is the total relative velocity between the water 
and the foil in the body Z direction, divided by the total relative velocity 


between the water and foil in the body X direction. 
X= ARCTAN “Ye (1-28) 


If the foil is not aligned with the body X axis, then the foil has a fixed 


angle of attack and the total angle of attack becomes: 


A total “Sgmgeg 1 (1-29) 


The side slip angle is shown in figure 1-5b and can be represented 


Dy: . 








= Ved 
on = ARCSIN - = = | (1-30) 
V US +Vej eo eh 
During actual craft Operations, the angles of and will be 


less than 10° permitting use of the small angle approximation. Also, the 
Vand W velocities can be considered very small in relation to the U 
velocity.. Applying these simplifications, the equations fore and G 


become: 


BY 





Body X-Axis 





Water Surface 


Body Z-Axis 





(a) 
ANGLE OF ATTACK 


Body’ X-Axis 


> 
— K, 
Body Y-Axis 
aa [Us +Vet We 


Vie 


Strut 


(b) 
ANGLE OF SIDE SLIP 


Facure j-5 


7a 





Wee 


xX, = U ate Ait ed (1-31) 
VV ‘ 
Bi = — (1-32) 
U 
D. EXPANSION OF FORCE AND MOMENT TERMS 


With the aid Of figure 1-4 and Table 1, the forces and moments for 
the craft are easily written. It is important to note, that when a numerical 
value is substituted for one of the alphabetical a the sign must 
be positive or negative depending on the relative position from the CG. 

One o% the assumptions made upon cOmmencing simulation was that 
the struts contributed zero lift and that the foils contributed zero side force. 
Therefore, all body Y motion originates at the struts and all body Z motion 
Originates at the foils. 


The force and moment equations become: 


Fy= Free tFeee t Face tFxme + Fyes + Fx ps t Fre tMgGy +Tx (1-33) 
Fy =Fyco+Fypg + Fysg + gy (1-34) 
Pz = Facet Fape +Face +t Fime +92 (1-35) 
L = (Fa ly)pe + (Fly) se (FyLz) ps ~( Fy L2)s5 “(Fy baves (1-36) 


M =~(F; Lx)oe ~ (Fahy) pe ~(Falx)se- (Fa lsyme + (Fx Laer 
+E La) pp tC Fhe) set(Fylayme t CFxbez)es 


+ (FL 2) pe + (Fy Lass t Ty Lar | (1-37) 


N= (Fytes t (Fylx)ps + (Fy Ly) 55 (1-38) 


ZL 





All of these forces have a hydrodynamic origin except for the 
thrust and gravity terms. The thrust is always associated with the body 
axis and requires no further expansion. The gravity terms can be simpli- 
fied using the small angle approximation for the Euler angles. (See 

summary Of equations A-31, 32, 33.) 
The general expression for hydrodynamic forces in water axes 


coordinates is 
Fé = Ge AC; = (1-39) 


where L, D, S represent lift, drag and side force respectively. The lift, 
drag, and side force coefficients vary as functions of angle of attack, 
angle of side slip, submergence, velocity, flap and rudder nosigiens. 

Ideally, mathematical expressions would have been developed to 
correctly depict the interrelation of all the variables which affect the 
force coefficients. However, hydrofoil technology has not advanced to 
a position which would yield such expressions. Some of the major prob- 
lems encountered in deriving a mathematical expression are: 

1) The Occurrence Of cavitation and ventilation which are 


completely unpredictable. 


2) Even more unpredictable is the cessation of cavitation and 
ventilation. 
3) The nonlinearity of the hydrodynamic coefficients. 


4) Lack of sufficient test data to completely describe all of 


the above. 
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II. HYDROFOIL SIMULATION 


Given the equations of motion, the next step was to formulate a 
program which would simulate the model. The computer language chosen 
was Digital Simulation Language or DSL. DSL was chosen over FORTRAN 
because Of its adaptability to control problems and the ease with which 
the user can insert engineering relationships with a minimum amount of 
programming. 

This computer simulation is designed for the HIGH POINT PC(H)-1 
which has a canard foil-strut configuration. Drawings of the craft and 
the foils ae shown in figures 2-1 and 2-2. The after foil is one unit, but 
in the equations, it is considered to be divided into three separate segments, 


i.e. mid panel and port/starboard outboard panels. 


The following assumptions were made prior to beginning simulation: 


lee Craft equations of motion are valid only for the foilborne 
mode. | 

De Weight of the craft remains constant. 

oa Hydrodynamic coefficients are based on fully wetted surfaces, 


i.e. no cavitation or ventilation. 


4. No constraint was placed on crafts heading. 

Se Foils are considered flat surfaces vice dihedral/anhedral. 
6. The craft has no automatic control system. 

ee The craft may be operated in calm water or in a seaway. 


Before the main program could be assembled, subroutines had to be 
written to find the hydrodynamic coefficients of foil lift and drag, and strut 


drag and side slip. Using the simulation curves in Ref. 3 as a guide, the 
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curves were stored as data points in the SETUP routine. The subroutines 
INTERP and INTRP1 are interpolation routines to obtain the proper coeffi- 
cient from the curves. For INTERP, values of angle of attack and flap 
angle are used to enter the curves and Obtain a value of foil drag. Angle 
of side slip and submergence are used to enter INTRP1 to obtain strut drag. 
The curve for the lift coefficient of a foil in the fully wetted region is a 
straight line so this coefficient is found by merely solving the equation 

for the straight line. Curves plotted from the subetee are shown in 
Appendix B. Listings of the subroutines are shown in the Computer Program 
section. The main program was then assembled using figure 2-3 as a refer- 
ence for data flow. A complete listing of the main program is shown in 


the Computer Program section. 


A. SIMULATION OBJECTIVES 

The Object of the simulation was to commence the runs with the 
craft foilborne and in a steady state condition. Once this condition was 
achieved, the step-response of the craft would be studied by applying 
step functions to excite motion along and about each body axis. That is, 
perturbations should be introduced to separately excite: 


. Motion along the X-axis 
; Motion about the X-axis 


: Motion along the Y-axis 


1 
Z 
2 
4. Motion about the Y-axis 
o Motion along the Z-axis 
6 


Motion about the Z-axis 
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After the step response was Obtained, a sinusoidal sea would be 


inserted to observe the crafts motion in a seaway. 


B. PROBLEMS ENCOUNTERED 
There were three majOr problems encountered during the simulation 
that are worthy of discussion: (1) the unknown stability of the craft, 
(2) specification of units used in the derivation of the equations, and 
(3) numerical quantities for thrust versus speed. 
1. Craft Stability 

One of the first assumptions made was that the craft would 
be stable in the foilborne mode with no control system. Starting from this 
assumption, the first consideration became that of balancing all of the 
forces and moments. A speed of 36 knots was chosen to insure that cavi- 
tation did not exist at any of the wetted surfaces and that the control sur- 
faces were zero. The fixed angle of attack was assumed zero. 

It was readily apparent that zero fixed angle of attack was 
an invalid assumption because the angle of attack was to be the control- 
ling parameter in stabilizing the model. Reference 5 listed the fixed 
angles of attack as 0 .332 forward and 0.0558 aft. These values were 
used, but the model became unstable at approximately 4.9 seconds. Only 
slight variations of angle of attack above and below a critical value caused 
the model to either raise out of the water or to ee This pointed out the 
strong dependence of both lift and drag upon the angle of attack and the 


problem of balancing the force and moment equations. Referring back to 
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equations 1-33 through 1-38, One can appreciate the complexity involved 
in balancing these equations as they are all interrelated. In order to 
simplify the problem at all, only the forces in the body X and Z directions 
and the moment about the Y axis were considered. Another simplification 
was that of considering the entire forward foil as a control surface. The 
Angle of eee: program was then written to solve for the angle of attack 
which minimized the force in the body Z direction and the moment about 
the body Y axis. Once the force and moment were minimized, the program 
then checked the force in the body X direction and adjusted thrust to mini- 
mize it. A listing of this program, which is straight FORTRAN, is shown 
in the Computer Program section. The angles of attack determined from 
this program were then inserted into the main program. Significant results 
were Obtained, however, the model remained unstable. This indicated 
that there was no feedback to dampen the error, so the model was function- 
ing as an Open loop system. Considering the output quantities to be linear 
accelerations and velocities, angular rates and angular positions, while 
the input quantity me reference height above calm water and Xo is the 
ee cured quantity which is fed back, the block diagram can be represented 
as shown in figure 2-3. No graphs are shown of this mode because of the 
instability of the model. 

Stability of the model without a control system was questioned 
at this point and it was determined that assumption #6, i.e. the craft has 
no automatic control system, was invalid. Reference 5 confirmed the fact 


that the model was only stable in three degrees of freedom without a control 
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system. Additionally, it was determined that when the model was con- 
strained to the pitch-heave-surge mode of operation it was still only 
marginally stable without a control system, At this point, a control sys- 
tem was added and it will be discussed in detail in Chapter III. 

Ze Specification of Units 

Considerable confusion resulted from the interpretation of 
the various references cOncerning the units used in defining Euler angles 
and angular rotations. This problem has been alleviated here by the 
specific definitions in the Table of Symbols. 

on Thrust vs Speed 

Since nO graphical or numerical data was available to show 
the thrust vs speed relationship, the computer was utilized to determine 
a value of thrust for a given speed. Once the Angle of Attack program had 
found the proper angle of attack to minimize the force and moment, the pro- 
gram then checked the force in the body X direction and adjusted the thrust 
to minimize this force. This produced an initial condition for thrust which 
was then inserted fen the main program. The main simulation program incor- 
erated a thrust trim routine which is shown on the main listing. 

The thrust trim routine compares the magnitude of the force 
in the X direction with the thrust and then tests the difference with an 
arbitrarily small fixed error. If the difference is greater than the error, 
then thrust is set equal to the value necessary to make the force in the X 


direction equal to zero. At the same time, the acceleration in the X 
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direction is also set equal to zero. If the difference is less than 
the error, then the thrust is held constant and the acceleration in the X 


direction is allowed to vary. 
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III. AUTOMATIC CONTROL SYSTEM 


Once it became apparent that the model would remain unstable 
without a cOntrol system, the Only alternative was to add the control sys- 
tem. As was previously stated, the problem will be greatly simplified by 
constraining the model to the pitch-heave-surge mode Of operation. In 
so doing, the control system can be limited to the forward foil. It was 
assumed that once the control system for the forward foil was operational, 


it could be easily extended to the after foil. 


A. HEIGHT SENSOR 

For this control system, a height sensor was incorporated which 
sensed the change in height of the CG above calm water. This change 
was utilized as a positional feedback signal to control the flap angle on 
the forward foil. This was accomplished in the simulation by defining the 
reference height of the CG for the particular speed and subtracting the 
actual height to obtain a difference called DELH. This difference was 


multipled by a gain factor to produce the variable FLAP defined as: 
FLAP = DELH*COEFH (3-1) 


Figure 3-1 shows the contents of the automatic control block in the pro- 
gram listing. Figure 3-2 is the block diagram of this system showing 
submergence height as the measured quantity being fed back to the 


cOntroller. 
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The amount of flap is governed by the unbalance of the equations 
of motion. Since the flap angle is directly proportional to both lift and 
drag it makes a significant contribution to the balancing of the equations 
of motion. A parameter study was made and determined that the optimum 
value for COEFH was 0.85. Figures 3-3, 4, 5 show the response of the 
model to a positive 5 deg/sec disturbance applied to Q for 0.2 seconds. 
It must be pointed out that with the configuration used, the flap was 
capable of moving instantaneously because an inertia term had not yet 
been added. 

With the model stable using position feedback, the next refinement 
was that of utilizing acceleration, velocity and position signals for feed- 


back. It was realized very quickly that this was an unrealistic approach. 


B. RATE FEEDBACK 

With re model stable using position feedback, the next refinement 
was that of utilizing acceleration, velocity, and position signals for feed- 
back. It was realized very quickly that this was an unrealistic approach 
in that two of the three feedback signals were measured at the CG, 
therefore, neglecting any pitch component. The equation for flap then 
becomes: 

FLAP = GAIN] DELH-GAIN3 W-GAIN4 W (g=2) 

where GAIN1, GAIN3, and GAIN4 are gain Beneemee: Figure 3-6 shows 
the contents of the automatic control] block for this system and figure 3-7 


is the block diagram. The model became unstable for all values except 
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Figure 3-1. Automatic Control Block for Height 
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zero, of GAIN3 and GAIN4. No graphs were plotted because of the 


unstability of the system. 


oF RATE FEEDBACK WITH ACCELEROMETER 

This system was implemented by using an accelerOmeter to sense 
motion at the forward foil. The accelerometer corrects the body vertical 
acceleration by the pitch rate translated to the forward foil. This pro- 


duces the true acceleration at the forward foil. The equation is given by: 
hoo = (FL/m)-Ly 4 pQ (3-3) 


Figure 3-8 shows the contents of the automatic control block for this 
system and figure 3-9 shows the block diagram. This control system pro- 
vided stable operation as long as GAIN3 was less than 0.02 and ae 
was less than 0.2. Figures 3-10, 11, 12 are graphs of the system 
response to the same disturbance stated previously. 

All three of the systems listed so far were not physically realizable 
in that the flap was provided instantaneous response to any error signal. 


This is not the case in a real system because the flap possesses a finite 


amount of inertia which produces a time lag in the response. 


D. RATE FEEDBACK WITH REAL POLE 

The instantaneOus response problem was taken care of in this system 
by adding an inertia term. This entails adding a pole to the root locus 
and placing it on the real axis in the left half plane. The addition of the 
pole helped to stabilize the system, but also increased the overshoot and 


and time lag. The time lag associated with the real pole caused the 
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Figure 3~ 8. Automatic Control Svstem With an 
Accelerometer 
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motion of the flap to lag the motion of the model and thus cause the flap 
to continue moving after the error had ceased. It was also noted that 
when the acceleration and pitch angle gain apg became too large, 
the system went unstable. Comparison of figures 3-13, 14, 15, with 
3-10, 11, 12, confirms these statements. 

The contents of the automatic control block are shown in figure 
3-16. Inthe block diagram of figure 3-17, the real pole term is incor- 


porated into the controller block. 


Er RATE FEEDBACK WITH COMPENSATOR 

In Order to minimize the effect of the inertia term, a compensator 
had to be included in the system. The inertia term caused a time lag in 
the response, therefore, a lead compensator had to be added to make the 
overall response approximate that of the instantaneous response. The 
term that = actually inserted was a computer artificiality because it 
took on the role of both real pole and compensator. The compensator of 


the simulation takes on the form (P S+1)/(P,S+1). In the root locus, 


] 
another pole was added, however, it is so far to the left that it was 
ignored. 

Comparison of figures 3-18, 19, 20 with figures 3-10, 11, 12 
shows that the compensator did reduce the overshoot and time lag assoOci- 
ated with the real pole. Close observation shows frat the response of the 


model with the lead compensator closely approaches the response of the 


model with an instantaneous flap. 
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The contents of the automatic control block are as shown in the 
main simulation program in the computer program section. The block dia- 


gram for the system is shown in figure 3-21. 
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III. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK 


An invalid beginning assumption, i.e. the craft is stable without 
a cOntrol system, caused the loss of a great deal of time early in the 
research. This was not a complete loss, however, as it necessitated 
the writing of the ANGLE OF ATTACK program to determine the angles of 
attack and thrust which balanced the equations of. motion. 

The data obtained from the ANGLE OF ATTACK program illustrated 
the tremendous effect slight changes in angle of attack make on the 
motion of the model. Since these aoe are so slight and occur in such 
short time intervals, it would not be possible for manual control. This 
necessitates an automatic control system that will provide the necessary 
Outputs to operate the model according to the given inputs. 

A stable model was obtained as a result of this research, however, 
it is somewhat limited. Using the present model and expanding it to 
include the after foils, the model would be completely stable in six degrees 
of freedom. 

Recommendations for further study include; 


P Incorporate a complete automatic control system which would 
include all three foils. 

2% Produce a simulation for the transition period between hull- 
borne and foilborne conditions. 

oie A study of the step response of the model in all modes. 

4, A study of the stability of the model in both regular and random 


seas. 
Sie A study of the turning characteristics. 
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fransformation of Velocities from Body to Barth axes 
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Hositions in Narthy Axes 
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melative Velocity Components 
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and side slip are calculated in radians. 
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Hydrodynamic Forces 1neiieie ees 
(Hach foil equation and strut equation must be repeated 
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MOLL Lift 
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Foil Drag 7 
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Strut Drag 
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